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^ Abstract: We investigate the differences between the NLO evolution 

kernels in the Curci-Furmanski-Petronzio (CFP) and Monte Carlo (MC) 
factorization schemes for the non-singlet case. We show the origin of these 
differences and present them explicitly. We examine the influence of the 
' -H ' choice of the factorization scale in the MC scheme (given by the upper 

phase space limit) on the evolution kernels in this scheme. 
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1. Introduction 

Results presented in this contribution are needed for the exclusive mod- 
elling of the next-to-leading-order (NLO) DGLAP jjj evolution of the parton 
distributions in the Monte Carlo, see refs. [21 O H] for general scope of the 
project. In this approach the NLO evolution of the parton distributions is 
done by the Monte Carlo (MC) itself, with the help of newly defined ex- 
clusive evolution kernels. The important point is that the commonly used 
MS factorization scheme is not well suited for defining the exclusive kernels 
for the above MC modelling. For this purpose one needs to define a new 
factorization scheme, referred to in the following as the Monte Carlo (MC) 
factorization scheme, see refs. [3l [2]. The calculations of the non-singlet 
NLO kernels (inclusive and exclusive) in this scheme have been presented 
in ref. [5]. The properties of the new exclusive evolution kernels were also 
investigated in refs. [H [7] . 

Why do we need a new factorization scheme? Traditional factorization 
theorems formulated in the early 1980's, see e.g. refs. f8| l9| TO] separate 
collinear singularities only after the phase space integration, which is quite 
convenient for the analytical calculations. In the MC approach, however, it 
is unacceptable. 

In defining the new MC factorization scheme classic works of Ellis, 
Georgi, Machacek, Politzer, Ross (EGMPR) [8j and Curci, Furmanski, 
Petronzio (CFP) [9] are used as a reference and a starting point. The 
definition of the new MC factorization scheme is not yet consolidated - it 
is still defined case by case, for more details see refs. O [5l [2j^ 

In the following we shall concentrate on the standard inclusive evolu- 
tion kernels in the MC factorization scheme, focussing on the non-singlet 
DGLAP evolution. For more than two decades of the collider experiments 
the "industry standard" in the data analysis has been the MS factorization 
scheme. It is particularly suited for calculating the hard process matrix 
elements and PDFs and hence widely used. One of the essential properties 
of the MS is that it uses dimensional regularization, very convenient in the 
analytical calculations but rather unfriendly for the MC simulations. It is 
therefore necessary to compare the inclusive NLO DGLAP evolution kernels 
(depending only on the light-cone variable x and flavor indices) calculated 
in the AIS and MC factorization schemes. This issue will be discussed in 
the following. 

In addition we shall show how the choice of the type of the factorization 
scale in the MC scheme influences the inclusive kernels in this scheme. In 
the MC scheme, contrary to the MS, where factorization scale n is a formal 



^ In the complementary approach of refs. [Ill 1121 113j soft singularities are resummed 
first and collinear resummation is added next. 
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parameter, factorization scale is given by a kinematical variable enclosing 
the phase space from the above. The so called evolution time variable used 
in the MC simulation is just a logarithm of this factorization scale. 

Let us concentrate on a subset of cut-diagrams with two real partons 
contributing to the non-singlet NLO evolution kernel, in any scheme using 
axial gauge, see Fig. [TJ 



■j\j'\r\J'\ry-\j-\r\j-\j^ 



(a) 



(b) 



(c) 

Fig. 1. A subset of Feynman cut-diagrams contributing to the non-singlet NLO 
DGLAP kernel: (a) ladder interference diagram, (b) double ladder diagram with 
counterterm, (c) gluon pair production diagrams: real and virtual. 



2. Notation and definitions 

Our main intrest in this contribution are 2-real emission cut-diagrams, 
see Fig. [T] For the parametrization of the four-momentum of two emitted 
partons we use Sudakov variables: 

fcf = aiP^ + ftn^ + A:f^, i = 1, 2, (1) 

where p is the four-momentum of the incoming quark and n is a light- 
cone vector (we use light-like axial gauge where = 0). In the massless 

hmit of QCD we have = kf = 0, which fixes A = -2^) = 2^^- 
Additionally we use symbol q for the off-shell momentum q = p — ki — k2 
and k for the sum k = ki + k2- 

Our preferred variables are related to polar angles of the emitted partons: 
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We refer to their modulus as the angular scale variables. They are related 
to rapidity via r^j = In |aj| + const. 

The inclusive kernel in the CFP scheme [9\(MS) is defined as twice the 
residue in e = of the bare PDF denoted as Id 



= 2Reso T, 

r = xPP 



(3) 



is a two particle phase space, 5 (^x — = 5(1 — x — ai — 02) defines the 

Bjorken x variable, C is a color factor, g is the strong coupling, W{ki, k2, e) 
originates from the 7-trace of individual diagrams contributing to F, is a 
formal energy scale of the MS scheme - the factorization scale, term /i^^ 
originates from dimensional regularization in n = 4+2e dimensions. Finally, 
the theta function @{s{ki, /C2) < Q) defines the upper phase space limit for 
two emitted partons using a function of kinematical variables s{ki, ^2). 

The definition of the inclusive kernel in the MC factorization scheme is 
similar: 

G = j d^5{x-'^ xCg^ W{ki,k2,e = 0) Q{Q > s{ki,k2) > go)- 

(4) 

In this case the upper phase space limit Q plays the role of the factorization 
scale. As we shall see, changing the type of variable used as the factorization 
scale implies modification of the MC scheme. Typically, the factorization 
scale (upper limit) is defined in terms of angular scale variable s{ki, /C2) = 
max{|ai|, |a2|} or transverse momentum s(fci, /C2) = max{|ki_L|, |k2_L|}. We 
shall refer to these two choices as angular ordering (a-ordering) and trans- 
verse momentum ordering (fe_|_-ordering). 

We also introduced a lower cutoff qo for the purpose of regularization 
of the overall scale singularity (in CFP it is done in a dimensional man- 
ner). The distributions W of the CFP and MC schemes coincide in four 
dimensions, W{ki,k2,e = 0) = W{ki,k2)- 

In the analytical integration over the phase space in the calculation of the 
inclusive kernels the parametrization of the phase space is in practice fully 
determined by the choice of the phase space upper limit. For instance, if the 



^ We skip the flavor indices as we are considering the non-singlet kernels, where nearly 
all diagrams contribute to the qq kernel. 
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angular variable is chosen, we also use it for the phase space parametrization: 

In the above an overall azimuthal angle dependence is already integrated 
out (only relative angle between emitted partons is kept). The same formula 
holds for the MC scheme in four dimensions (e = 0). In case of A;_L-ordering, 
the analogous parametrization reads: 



,,T, 1 ^l+2t daida2^ , ,l+2e,l+2e (a\ 

"^^'^ = 4(2^0^^^ ^2! • (6) 



3. Single-logarithmic divergent diagrams 

In this section we show the general structure of single logarithmic diver- 
gent diagrams contributing to the NLO evolution kernel using the example 



of the ladder interference diagram of Fig. 1(a), which we shall call Bx in 
short. We analyze the integration procedure in the MS scheme and in two 
MC schemes (with factorization scale in terms of angular scale variable and 
transverse momentum). 

Choosing the upper phase space limit in terms of angular variable, the 
bare PDF (kernel) contribution of Bx diagram in the MS factorization 
scheme reads: 

X / (0:102) 0(1 — X — 01—02) / dih+2e 

J oi 02 J 

X / daida2-^, -rB^.(ai/a2, 6*, oi, 02, e) G( max{ai, 02} < Q) >, 

Jo r (01,02) 'J 

(7) 

where Tbx is a dimensionless function and cp' = a^ + ^-^^a^ + 2ai • a2. 
As it was already pointed out, Bx diagram features only single collinear 
divergence (single ^ pole). This singularity can be isolated in the very be- 
ginning of the calculation (integration over the overall scale). It is done 
by means of introducing new integration (scale) variable using the iden- 
tity Q{Q > s{ki,k2)) = jiQ dQ5{Q = s(A;i,A;2)), using at the same time 
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dimensionless variables = ^ Eq. M reads: 

"^"^ '^"^ iaia2f'Sil - X - ai - as) ^ dQQ^'~^ J d^i+2e (8) 



X / dyidy2^^^^P^ — -TBx{yi/y2,6,Oii,a2,e) 6{l - max{yi,y2})\. 
Jo r (2/1,^2) 'J 

The overall scale singularity is now factorized in the form of the integral 
/o^ dQQ^^~^ = If we perform the same calculation in the 4 dimensional 
MC factorization scheme, (e = 0) the scale singularity is regularized by 
means of a cutoff 0(max{ai,a2} > qo). Then the integral containing the 

scale singularity takes the following form J^^ ^ = ln(Q/go)- 

Since Bx diagram does not features additional coUinear singularities we 
can perform the action of the pole part operator of the CFP scheme imme- 
diately and evaluate the remaining integrals in four dimensions. However, 
there are still infra-red (IR) singularities originating from the integration 
over light-cone variables a. In the CFP scheme they are regularized in 
geometrical (non-dimensional) manneij^ and they do not contribute any ad- 
ditional collinear poles 

X / dyidy2 ^,,^ ^ TBx{yi/y2,d, ai, a2, e = 0) 6(l - max{yi,y2}) ■ 

Jo r{yi,y2) 

(9) 

In the MC scheme gets replaced by ln((5/(7o)- It means that the inclusive 
(integrated) kernel contributions from Bx diagram (and also from all dia- 
grams with only single collinear divergence) are the same in both the MS 
and the MC scheme with angular factorization scale. A slight difference in 
the overall normalization is a matter of convention: 

P{x) = 2T(x). (10) 

What will happen if we change the type of the upper phase space limit? 
In ref. 171 it was checked for the MS scheme that the inclusive NLO kernels 



The IR singularities in the CFP and MC schemes are regularized by principal value 
prescription: ^ — >■ ■ We also use the following notation of CFP for divergent 

integrals: da^^^ = Jo and dalna^^^ = h. 
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in this scheme do not depend on the choice of the upper hmit (as expected). 
The question is now whether the same holds in the MC factorization scheme. 

In the case of single logarithmic divergent (single ^ pole) diagrams, like 
Bx, it is true again. It can be seen without performing the actual cal- 
culation - just by means of inspecting once again the calculation in the 
MS scheme. We could write equation analogical to Eq. ([T]), but with the 
upper limit in terms of another variables, for instance using transverse mo- 
mentum G( max{/ci_|_, A;2_l} < Q). Next we would introduce dimensionless 
variables y'^ = ^ and isolate once again the overall scale singularity. Since 

we consider diagrams free from internal ^ poles we can take the residue or 
differentiate over ln(Q), obtaining the same result for both the MS and MC 
schemes, independently of the type of the variable used in the MC scheme. 

Summarizing, we have shown that the contribution to the NLO evolution 
kernel from the diagrams with single logarithmic divergence gives the same 
kernel contributions in all considered factorization schemes: the MS, MC 
with angular and transverse momentum factorization scale. Moreover the 
same argumentation will be applicable to the single logarithmic parts of the 
other diagrams. However, their double logarithmic parts will deserve special 
considerations, see below. 



4. Double-logarithmic divergent diagrams 

In this section we consider the inclusive kernel contributions of the 
double logarithmic divergent (double pole) diagrams in the MS and MC 
schemes. We use the example of the double ladder diagram (called Br) and 



and virtual) of Fig. |l(c) (called Vg 



its counterterm of Fig. 1(b) and the gluon pair production diagrams (real 



The double divergent diagrams are the sources of differences between 
inclusive kernels of MS and MC schemes. The reason is that MS kernels 
contain additional terms, originating from the artefacts of the dimensional 
regularization, generally form mixing e x ^ = ^ coming from additional e 
terms in the phase space (kinematic) or in the 7-trace (spin). 



4-1. Double ladder diagram and its counterterm 

The structure of the kernel contribution of the double ladder graph (Br) 
is analogical to the one of the ladder interference diagram given in Eq. 0. 
In case of the counterterm (Ct) the situation is slightly different due to 
the presence of additional projection operator. The calculations in case of 
Br and its counterterm are more complicated and we do not present them 
explicitly in this contribution. We will present only the final results. For 
the details of the calculations we refer the reader to ref. [5] . 
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The inclusive kernel contribution of the double ladder graph and its 
counterterm in the MS factorization scheme reads: 



'yBr-Ct{x)=[^)^Cl 



1 



■ X 



+ 3(1 



1 + 



ln(a;)) 



- 4/o - 41n(l - x) + 2\n^{x) 



{l + x)[ \n{x) + ^ln^(a;) 



:iii 



It does not depend on the choice of the upper phase space hmit. The same 
kernel contribution in the MC scheme with angular factorization scale is the 
following: 



l + x^ 
1-x 



87o-81n(l-a;) + 41n2(x 
+ {l-x)(6-2 ln(x)) - (1 + x) (2 ln(x) - ln2(x)) 



(12) 



and in the MC scheme with transverse momentum factorization scale: 



^Br-Ct 



1 +X^ 



-8I0 -81n(l -x) +41n2(x) 



+ {l-x)(6 + 2 ln(x)) - (1 + x) (2 ln(x) + ln2(x)) | . 



(13) 



The difference between the MS and MC schemes originates from the prod- 
uct of double pole ^ and the coefficient terms expanded to 0{e^). There 
are two kinds of these terms. The first are from the phase space and depend 
on the type of the upper phase space limit (factorization scale in the MC 
scheme). The second are from the 0{e^) terms in the 7-trace and do not 
depend on the type of the factorization scale in the MC (spin part). 

The difference between the MS and MC kernels written explicitly reads: 

P|,_Ci(x) - 2TBr-Ct{x) = 7t_cM) + nT-Ct{x), (14) 

where index s indicates the type of the factorization scale in the MC scheme. 
The spin part of the difference is the same for all choices of factorization 
scales in the MC and reads: 

(15) 



g) C|,4(l-x)ln(x). 



The kinematical difference is due to the choice of the upper phase space 
limit in the MC scheme. For the angular factorization scale it reads: 



V27r/ 



2(1 + x) ln2(x) - 4(1 - x) ln(x) 



(16) 
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For the transverse momentum factorization scale it vanishes: 



:i7) 



The kinematical difference is due to additional phase space factors in n- 
dimensions. For angular phase space of Eq. ([s]) it is the term (aia2)^'^, which 
after expansion, multiplied by the double pole, contributes 2eln(aia2) x ^ 
to the NLO DGLAP kernel. In case of the transverse momentum phase 
space of Eq. ^ this kind of term is missing, i.e. (^) ~ 

The kinematical difference J'^r-ct('^) ^^^^ exactly equal to the dif- 
ference between kernel contributions in two MC schemes, one with angular 
and another with transverse momentum factorization scales: 



PBr-Ct{x)~PBr-Ct(^) — (2^) 



2(l+rc)ln2(x)-4(l-x)ln(x) 



(18) 



The above indicates that the MS factorization scheme is affine with the 
MC scheme with transverse momentum factorization scale. Maximum k± 
is effectively representing a formal scale parameter fi (factorization scale) 
of the dimensional regularization in the MS. As was shown in ref. [7] MS 
kernels do not depend on the type of the upper phase space limit. If we 
choose phase space enclosing in form of another kinematical variable (not 
k±), the additional ~ e x ^ terms in the MS intervene and self-correct 
back to the case of transverse momentum. In case of the 4-dimensional 
MC schemes there is no automatic self-correcting mechanism like in the 
MS and kernels in the MC schemes with different types of the upper phase 
space limits will differ. 

Summarizing, the "spin difference" between kernels in the MS and MC 
schemes is the same, whereas the "kinematical difference" depends on the 
type of the upper phase space limit (factorization scale) in the MC scheme. 
For instance, for the factorization scale in terms of the light-cone minus 
variables |kj_L|/y^ the kinematical difference is due to the term (0102)"^ 
and reads: 



Kin |kix|/ 
Br-Ct 



;i-h2;)ln2(x)-2(l-x)ln(x) 



(19) 



The case of virtuality {—q^) upper limit was analyzed as well. The 
general result is the same, but the details of the calculations are a little bit 
more complicated. 



4-2. Gluon pair production diagrams 
The 2-real gluon pair production diagram (Vg) of Fig. 



l(c)| feat 



ures 



internal collinear singularity due to the production of a pair of collinear 
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gluons. Analogous diagram with production of a quark antiquark pair also 
contributes. As has been shown in ref. IlTJ, also in this case the MS kernels 
do not depend on the choice of the upper limit. It is ensured (self-correction 



mechanism) by the corresponding virtual diagram of Fig. l(c^ 



The question is what happens in the 4-dimensional MC scheme? The 
internal collinear singularity leads to the double e pole, which suggests that 
there could be a difference between kernel contributions in the MC schemes 
with different factorization scales. On the other hand, there is no additional 
projection operator as in the case of the double ladder graph and its coun- 
terterm, which suggests otherwise. In this case it is the virtual graph, which 
ensures the independence of the MS kernel contribution on the type of the 
upper limit. 

The real and virtual graphs are combined - this combination in the MS 
and MC schemes is done in the same manner. Partial cancellation of the 
^ pole occurs, with the remnant giving rise to the running of the coupling 
constant. It originates from the gluon self energy graph (the virtual graph 
in Fig. 1(c)) from the integration of the term like ln(Q//xij), where is 
the renormalization scale and Q is a kinematical variable to be chosen - 
the argument of the running coupling. The different choices of the running 
coupling argument Q do not affect the MS kernels, but they may affect the 
kernels in the MC schemes. The results of the preliminary studies indicate 
that if we choose the variable Q in terms of transverse momentum both the 
MC and MS schemes give the same results. 



5. Conclusions 

We investigated the differences between the non-singlet NLO DGLAP 
evolution kernels in the Curci-Furmanski-Petronzio scheme {MS) and in the 
factorization scheme defined for the purpose of MC simulations. We show 
that the origin of these differences are due to the different definition of the 
projection operator in both schemes and examine explicitly the source of 
the differences. 

We also examined the differences between MC schemes with factoriza- 
tion scales given in terms of rapidity related variables and transverse mo- 
mentum, indicating the source of the differences explicitly. The type of the 
factorization scale in the MC scheme given by the upper phase space limit 
is related to the evolution time variable in the MC simulations. Knowing 
the differences between the kernels for different factorization scales is criti- 
cal for understanding basic features of the MC factorization scheme and in 
particular what happens when switching between different evolution time 
variables in the MC simulation. 
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